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(54) iVIethod and apparatus for classifying heartbeats in an ECG waveform 



(57) A method and apparatus for classilying heart- 
beats obtains ECG data from a plurality of ECG wave- 
forms vvhicn are in turn obtained from signals received 
from a plurality of ECG electrodes QRS detection logic 
delects heartbeats in the ECG data Classification logic 
classifies tfie delected heartbeats into categones based 



on shape and/or timing. This classification is done by 
comparing each heartbeat against a group of templates 
corresponding to one or more heartbeat classifications. 
The templates are updated to track changes in the mor- 
phology of the heartbeats. The heartbeat classification 
IS displayed to the cardiologist or medical professional 
(or diagnosis of the condition of the parent's heart. 
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Description 

Cross Reference to Related Applications 

5 This application is related tolhe following commonly assigned patent applications, filed on even date herewilii: 



Titie 


Serial Number 


Mettiod And Apparatus For Detecting Heartbeats In An ECG Waveform 




Method And Apparatus For Calculating A Heart Rate In An ECG Waveform 




Method And Apparatus For Creating A Representative Heartbeat From An ECG Waveform 




Method And Apparatus For Obtainmg Heartbeat Measurements From An ECG Waveform 





Field of tlie Invenllon 

This invention relates to the electronics c; -Iry field More particularly, this invention is a method and apparatus 
for classifying heartbeats in an ECG waveform. 

^0 Background of the Invention 

Designers of medical insirumeniation such as cardiographs face many difficult challenges in their jobs. The devices 
Ihey design are expected to deliver high quality information about the electrical activity of a patient's heart to a cardi- 
ologist or other medical professional, so that a correct diagnosis of the condition ol the patient's heart can be made. 

^"5 Unfortunaiety. the ECG electrodes connected to a patient usually deliver ECG data to a cardiograph that comprise not 
only information showing the electrical activity of the patient's heart, but also electrical noise This noise can make up 
most of the ECG data, and can corrupt and totally overwhelm the portion of the ECG data that contain information 
about the electrical activity of a patient's heart This problem is especially acute in hostile environments, such as a 
patient undergoing a stress or exercise test, where (he noise can be quite extreme Unless the medical instrumentation 

^0 designers are successful at designing a medical instrument that analyzes this ECG data to eliminate or reduce the 
effects of this noise the cardiologist or other medical professional will find It difficult, if not impossible, to obtain infor- 
mation about a patient's heart, such as the classification of Ihe type ol heartbeat, that is usable in arriving at a correct 
diagnosis of the condition of the patient's heart. 

Summary of the Invention 

A method and apparatus for classifying heartbeats obtains ECG data from a plurality of ECG waveforms, which 
are in turn obtained from signals received from a plurality ol ECG electrodes. QRS detection logic detects heartbeats 

•n the ECG data Classification logic classifies the detected heartbeats into categories based on shape and/or timing. 
■^0 This classiricalion is done by comparing each heartbeai against a group of templates corresponding to one or more 
heartbeat classifications The templates are updated to track changes in the morphology of the heartbeats, The heart- 
beat classification is displayed lo the cardiologist or medical professional for diagnosis of the condition ol the patient's 
heart 

-^s Description of the Drawings 

Fig. 1 shows a block diagram of the medical device of the preferred embodiment of the invention. 
Fig. 2 shows a block diagram of the medical device of Fig. 1 in more detail. 
Fig 3 shows the processor ol the computing unit of Ihe medical device of Fig. 2 in more detail. 
Figs 4 ana 1 2 show flowcharts of the operation of th-^ *ardiograph of the preferred embodiment of the invention. 
Fig. 5 shows a high level flowchart of the operatior the computing unit of the preferred embodiment ol the 
nvention 

Fig 5 shows a flowchart ol the operation of the GPS detection logic ol the preferred embodiment of the invention 
Fig 7 shows a flowchart of the operation of ;he heart rate calculation logic of the preferred embodiment of the 
^= invention 

Fig 5 shows a flowchart of the operation ol ihe classification logic of the preferred embodiment of the invention 
Fig 9 shows a flowchart of the operation of the alignment logic of the preferred embodiment of the invention. 
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Fig 10 s^cws a Howcnart of the operdUon of the representative heartbeat creation logic of the preferred embodi- 
f^ent of (he invention 

Fig 11 shows a flowchart of the operation of the measurements logic of the preferred embodiment of the invention . 
Fig 13 shows a graph of three exemplary ECG waveforms used by the QRS detection logic of the preferred 

ombodimenc of the invention 

Fig 1 4 shows a graph of an exemplary activity function used by the QRS detection logic of the preferred embod- 
irneni of the mventton 

Fig 15 shows a graoh of exemplary classified heartbeats. 

Fig IS snows a graph of exemplary heartbeats being aligned by the alignment logic of the preferred embodiment 

of the invention 

Fig 1 7 shows a graph of exemplary aligned heartbeats being timesliced by the representative heartbeat creation 
logic of the preferred embodiment of the invention. 

Fig. 15 snows a piintout or display of a representative hearlbeat without measurements. 
Fig. 1 9 shows a printout or display of a representative heartbeat with measurements. 

Detailed Description of the Preferred Embodiment 

I. Overview 

Fig. 1 shows a block diagram of the medical device of the preferred embodiment of the invention Medical device 
10 comprises acquisition unit 20. electrodes 25 cardiograph 40. and computing unit 60. In the preferred embodiment, 
cardiograph ^0 and acquisition unii 20 are separate components of a PageWriier XLi. manufactured by the Hewlett- 
-Packard company, modified to execute the flowcharts of Figs. 4 and 12 of the preferred embodiment of the invention 
Computing unit 60 is a HF Vecira personal computer, suitably programmed to execute the flowcharts of Figs. 5-11 of 
the preferred emt>ocl[ment of the invention 

Fig 2 shows a block diagram of medical device 1 0 in more detail. Cardiograph 40 contains acquisition unit interface 
41 processor 45 printer 47 and computing unit interface 49 Processor 45 executes the flowcharts of Figs. 4 and 12 
of the preferred embodimenl of the invention Computing unit 60 contains cardiograph interface 61. processor 65 
display 66. input device 67 memory 63 and storage 69 Processor 65 executes the flowcharts of Figs. 5-11 of the 
preferred embodiment of the invention While Fig. 2 shows medical device 10 as containing discrete components, 
mose skilled in the art will appreciate thai medical device 1 0 could be a single unit that contains each ol the components 
shown in Fig 2 or contain a different number of discrete components, and still fall within the spirit and scope of the 
invention 

Fig. 3 shows processor 65 of computing unit 50 of medical device 10 in more detail. Processor 65 contains QRS 
deieciion logic 71 heart rate calculation logic 73 classification logic 74. alignment logic 75. representative heartbeat 
creation logic 77 and measurement logic 75 In the preferred embodiment, each of these logic blocks is performed by 
software written to perform the functions of relevant portions of the flowcharts shown in Figs. 5-11. and this software 
;s executed by processor 65. Alternatively, some or all of logic blocks 71 -78 could be special purpose hardware, such 
as contained in an application specific integrated circuit, designed to perform functions of relevant portions of the 
flowcharts shown in Figs 5-11 

^iq. 4 shows a flowchart of the operation of cardiograph 40 of the preferred embodiment of the invention. In block 
101 ECG signals are received from electrodes 25 of acquisition unit 20. In the preferred embodimenl. these signals 
are digital signals sampled at a high sampling rale. Block 103 decimates and filters the sampled ECG signals outside 
cl a predeternnined bandwidth In the preferred embodiment, the predetermined bandwidth is .01 Hz to 150 Hz. and 
(he decimation process reduces the number ol samples to one eighth of the number of original samples. Block 105 
makes twelve ECG leads from the ten elecirodes m a conventional manner The signals contained on the twelve ECG 
leads will be referred to herein as "ECG waveforn^s* and the information contained thereon will be referred to herein 
as "ECG data" Those skilled in the art will appreciate that the actual number of electrodes or leads may be different 
than discussed aoove and still fall within the spirit and scope of the invention Block 108 sends the ECG data on the 
ECG waveforms io computing unit 60 The flowchan ends in block 109. 

"^ig 5 shows a high level flowchart of ihe operation of computing unit 60 of the preferred embodiment of the 
invention Block 201 receives Ihe ECG waveforms containing the ECG data from cardiograph 40. Block 203 forward 
fibers the ECG data In the preferred embodiment this filter is a high pass filter used as part of a forward/reverse 
filtering scheme to remove baseline wander while presen/ing low frequency information in the ECG data. Block 300 
calls a subroutine thai detects the heartbeats (i.e QRS complexes) in the ECG waveforms This subroutine calculates 
rin dcli\/iiy funcijon from a subset of the ECG waveforms determined to be least noisy and uses this activity function 
to search for t'-eartbeats This allows for true heartbeats to be detected while discarding false "noise" beats. The op- 
eration ol ims subroutine will be descnbed in more detail later m conjunction with the discussion ol Fig. 6. 
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Block 400 calls a subrouiine that calculates the patient's heart rate. This logic determines the inle/vals between 
the heartbeats discards a percentage oi the shortest and longest intervals, and averages the remaining inten/als to 
arrfVG at the patient's heart rale. This results m a robust calculation of the heart rate even m the presence of noise 
falsely delected as heartbeats and missed beats common m noisy environments. The operation of this subrouiine will 
be cescribed in more detail later in con)unciion with the discussion ol Fig. 7 

Block 500 calls a subroutine that classifies heartbeats. This classification is done by comparing each heartbeat 
against a group of templates corresponding to one or more heartbeat classificalions. The templates are updated to 
track changes in the morphology of the heartbeats. The operation of this subroutine will be described in more detail 
later in coniunclion with the discussion of Fig. 8. 

Block 205 reverse filters the ECG data. In the preferred embodiment, this filter is a high pass filler used as part of 
a forward/reverse filtering scheme lo remove baseline wander while preserving low frequency information In the ECG 
data 

Block 600 calls a subroutine that aligns heartbeats prior lo representative heartbeat creation. This logic slides the 
heartbeats across an alignment template heartbeat lo calculate when the heartbeats are aligned, and performs ad- 
justments to reduce the effects of noise or |itler on the different ECG waveforms. The operation of this subroutine will 
be described in more detail later in conjunction with the discussion of Fig. 9. 

Block 700 calls a subroutine that creates a representative heartbeat from the aligned heartbeats. This logic time 
slices through the aligned heartbeats, discarding a percentage of the smallest and largest magnitudes of the aligned 
heartbeats at each instance of lime and averaging the remaining magnitudes lo produce a representative heartbeat. 
This trimmed averaging technique results in a high quality represenlaiive beat, since samples from noise and misclas- 
sified beats are discarded, The operation of this subroutine will be described in more detail later in conjunction with 
the discussion of Fig. 10. 

Slock 600 calls a subroutine thai measures various aspects of a representative heartbeat. This logic analyzes the 
representative heartbeats from a group of ECG waveforms to determine an earliest QRS onset and latest QRS offset, 
and uses these values to perform a variety of measurements. This results in robust measurements even in very noisy 
environments The operaiion of this subroutine will be described in more detail later in conjunction with the discussion 
of Fig 11 

Block 210 displays the representative heartbeat created by subroutine 700 and. optionally the measurements 
obtained by subroutines 800 and 400. on display 66 of computing unit 60 (Fig. 2). Examples of these displays are 
shown in Figs 18 and 19 Block 220 sends the representalive heartbeat and measurements for each ECG waveform, 
including the heart rate measurement calculated by subrouiine 400. back to cardiograph 40, Cardiograph 40 processes 
this information in accordance with the flowchart of Fig, 12. The flowchart ends in block 249 

Fig 12 shows how cardiograph 40 processes the information received from compuling unit 60. Block 150 receives 
the representative heartbeat and measurements, including the heart rale measurement, sent by block 220 of Fig. 5 
Block 190 prints the representalive heartbeat created by subroutine 700 and. opUonally. ihe measurements obtained 
by subroutines 800 and 400. on printer 47 of cardiograph 40 (Fig. 2). Examples of these printouls are shown in Figs. 
IS and 19 



II. QRS Detection 

JO 

Fig 6 shows a flowchart of the operation of subroutine 300. performed by QRS detection logic 71 of compuling 
unit 60 of the preferred embodiment of the invention. Block 302 performs an initialization process that learns about the 
ECG data More specifically, the first time through the subroutine this process analyzes the first few seconds of ECG 
data to determine a preliminary interval between peaks and the magnitude of an average peak. During routine process- 
-5 tng i\.e subsequent limes through the subroutine^ block 302 continues to update the information obtained the first time 
through (he subroutine. Block 301 receives ECG data from three ECG waveforms out of the twelve ECG waveforms 
received by computing unit 50 in block 20i (Fig. 4) A graph of three exemplary ECG waveforms is shown in Fig. 13 
in the preferred embodiment block 301 selects these three ECG waveforms as the three ECG waveforms thai have 
optimal noise characteristics. This is done by continuously calculating the signal quality on the twelve ECG waveforms 
and ranking these waveforms from highest to lowest signal quatiiy 

Block 303 calculates an activity function out of the three leads of ECG data An activity function is a signal math- 
am;4tically derived from the ECG data which emphasizes characteristics of the heartbeat while minimizing the influence 
of noise to enaole more accurate heartbeat detection In the preferred embodiment, the activity function is created by 
decimating ihe ECG data by one half, filtering the data using a bandpass filter and taking the absolute value of the first 
difference The absolute first differences from the three ECG waveforms are then summed together, threshotded and 
smoothed lo improved noise performance. Noise statistics (i s signal quality) across all twelve ECG waveforms are 
calculated and updated A graph of an exemplary activity function is shown in Fig. 14 

Block 305 checks Ihe noise statistics to see if unaccepiably high noise was delected. If so. block 306 issues a 
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'high noise" message, which is displayed on display 66 of computing unit 60 {Fig 2). This error causes the subroutine 
to terminate abnormally 'oy ending m block 396. In the absence of high noise, block 305 is answered negatively. Block 
309 looks to see if it is to detect another heartbeat, if not ihe subroutine returns m block 399 to block 400 of Fig 5 
If block 309 is answered affirmatively, subroutine 300 pertcrms up to three different types ol searches to find each 
s heartbeat. The first search is an on-time search This search is normally used to detect heartbeats that appear wilhin 
a snnall window of Iheir expected time If the on-iime search fails to detect a heartbeat, a modified on-time search is 
used The modified on -time search is similar in operation to the on-time search, but can detect heartbeats in low noise 
environments that the on-time search might miss, such as early beats. If both the on-time search and the modified on- 
time search fail to detect a heartbeat, a late beat search is performed. This search delects heartbeats that appear later 
^0 than their expected lime 

The operation o( the on-time search will now be discussed in more detail. Block 310 performs an on-time search 
of ihe activity function for QRS complexes in the ECG data. In the preferred embodiment, block 310 searches over 
1 1 5°o of the current interval length between heartbeats (as determined in block 302) with a threshold that starts at 60% 
of the average peak value of ihe activity function The threshold is then linearly decreased over time to 40**'o at the end 
'5 of the search window. Any local maxima discovered above this linearly decreasing threshold is considered a "peak". 

Block 315 checks to see if any peaks were detected. If so block 320 selects the true heartbeat {i.e.. QRS peak} 
from the detecteo peaks by looking at the limes each peak occurred. Those skilled in the art will appreciate that the 
^bove-described search could detect multiple peaks in a noisy environment all but one of which would be noise The 
peak thai occurred closest to the time the next heartbeat was expected is considered a true heartbeat, and timing 
20 information about this heartbeat is stored in storage 69 ol computing unit 60 (Fig. 2). Flow of control loops back to 
block 309 to look for another beat to detect. 

Referring again to Fig. 6 the operation of the modified on-time search will now be discussed If block 315 is 
answered negatively block 330 checks to see if the noise level is low. If so. block 335 searches the activity function 
by using a modified on-iime search In the preferred embodiment this search is performed by using a constant lower 
threshold such as 54% of the average peak acltvtly (unclion value over 1 1 5% of Ihe current interval length. Block 335 
checks to see if any peaks were delected, if so block 333 selects the true heartbeat from the detected peaks by 
selecting the peak that occurred closest to the time the next heartbeat was expected. Information about the heartbeat 
selected in block 338. such as beat timing {nformation. is stored in storage 69 of computing unit 60 (Fig. 2). Flow of 
control loops back to block 309 to look for another beat to detect 
^0 The operation of the late beat search will now be discussed. If either blocks 330 or 335 are answered negatively 

block 350 performs a late beat search In the preferred embodiment, this block searches the activity function over three 
R'R pnlerval lengths using a linearly decreasing threshold, as was done with the on-time search. Block 355 checks to 
see pf any peaks were delected. If so block 35S selects the true heartbeat from the detected peaks by selecting the 
{■rs( peak that it finds Information about the heartbeat selected in block 358. such as beat timing mfomialion. is stored 
-•5 in storage 69 of computing unit 60 (Fig 2) Ffow of control loops back to block 309 to look for another beat lo detect. 
If block 355 IS answered negatively block 370 issues a "detector laiJed" message that is displayed on display 66 of 
com^puung unit 60 Since no beats were detected this causes an abnormal termination of the subroutine in block 396. 



III. Heart Rate Calculation 

Fig. 7 shows a flowchart ol the operation of subroutine 400. performed by the heart rale calculation logic 73 of 
computing unit 60 of the preferred embodiment of the invention Block 40i sets a counter called RR_ctr = 0. Block 403 
roads and stores the R-R interval between the hrsi heartbeat delected in the ECG waveform and ihe second heartbeat 
detected in the ECG waveform In the preferred embodiment this detection is done using information stored by ORS 
detection logic 7i in subroutine 300. although ccnvenlicnal methods of delecting R-R intervals in an ECG waveform 
could also be used Block 405 increments RR_ctr Block 410 checks lo see if the counter is less than a maximum 
counter value and there are additional heartbeats stilt available for analysis in the ECG waveform. If both of these 
conditions are true flow of control loops back tob»ock403 until one condition is no longer true. An alternate embodiment 
has been contemplated where block 403 is a timer and where block 410 checks to see if a maximum time has elapsed 
P"or examole. if a maximum time was set for i0 seconds only heartbeats occurring in the most recent 10 second period 
of :ime are used to calculate the heart rate 

When block 41 0 is eventually answered negatively block 420 checks to make sure at least a minimum number of 
heartbeats were analyzed by the looo made up of blocks 403-410 

If block 420 is answered affirmatively block 425 sorts the R-R intervals from shortest lo longest. Block 430 then 
discards a peicenlage of the shortest and longest -R-R intervals In a noisy environment, the QRS detector may incor- 
rectly detect noise as -h heartbeat and m,ay inccfroctly miss a real heartbeat. These errors result in incorrect R-R 
inten/als boih too short and too long. The trimmed ^-iverage done here results in a robust and accurate heart rate 
CHiCulation even in Ihe presence of false detections nod misseo oeats Wiih low noise and/or arrhythmias, this trimmed 
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average technique also produces an accurate calculation of the heart rate. 

I n the preferred embodiment, block 400 discards 2S% of the shortesi and 25% of the longest R-R intervals, although 

different values could be used. Block 435 then averages the remaining R-R intervals. Block 440 converts this average 

R-R interval to a heart rate. Block 445 smoothes the heart rale determined by block 440 by averaging it with a prede- 
5 lermined number of past heart rates. In the preferred embodiment, block 445 averages the current heart rate with the 

past two heart rates In any event, ihe heart rate determined by block 445 (or by block 440. if the smoothing step of 

block 445 is not desired) is stored by block 450 in storage 69 of computing unil 60 (Fig. 2). The subroutine returns in 

block 499 to block 500 of Fig. 5. 

If block 420 is answered negatively block 460 simply computes the average R-R interval of the small number of 
70 R-R intervals that were read and stored in block 403 This average R-R inten/al is converted to a heart rate in block 

440. and the heart rate is stored by block 450 in storage 69 of computing unit 60 (Fig. 2). As before, the subroutine 

returns in block 499 to block 500 of Fig. 5. 

IV. Heartbeat Classification 

Fig. 3 shows a flowcharl of the operation subroutine 500. performed by classification logic 74 of computing unit 
60 of the preferred embodiment of Ihe invention In the preferred embodiment, subroutine 500 is used to classify beats 
as "D " (dominant) "V" (ventricular ectopic) "S" (supraventricular ectopic), or "Q" (questionable), although other clas- 
sifications could be used 

'-0 Block 501 checks to see if there are more heartbeats to classify. If so. Block 502 gets beat timing information for 

the next heartbeat In the preferred embodiment, this information is obtained from QRS detection logk; 71 in a manner 
that has already been discussed, although conventional methods of obtaining this information could also be used. 
Block 503 normalizes the portion of the activity function (from QRS detection logic 71 or conventional means) around 
the beat to be used in the template matching process. Block 505 performs timing and physiologic measurements on 

?s both the activity function and on the beat to classify. These measurements are used in blocks 510 and 550 to assist 
in beat classification, as will be discussed later 

Block 510 checks to see if all measurements made by block 505 are wilhin physiologic limits (e.g. within a prede- 
termined width and height), if not. the beat is classified as questronable ("Q"). due to noise, and flow of control loops 
back to block 50i to see if there are more beats to classify. If block 510 is answered affirmatively, block 515 sets a 

'0 template counter to i . initializing a series of steps that compares the beat to classify with one or more templates. Block 
515 checks lo see if a counter called Num.Jemplales = 0. If so no templates have yet been created for this ECG 
waveform and flow of control skips down to blocks 539 and 540 to create a new template by saving this beat as the 
first template. Block 542 then uses timing and physiologic information to classify this beat and template as either "D" 
iDominant) "V" (ventricular ectopic) "S" (supraventricular ectopic) or "Q" (questionable). fVlost commonly this beat 

5 will be classified as 0 and the template classified as D/S lor "Dominant/supraventncular ectopic", since the vast majority 
of beats classified will be classified this way. and since both D and S beats have the same morphology and thus would 
match the same template, but vary by timing information ■- the S beat being earlier than the D beat. The method of 
atstinguishing between D and S beats is performed by blocks 545 lo 553. as will be discussed m more detail later 
Those skilled in (he art will appreciate that more than one template can be classified as "D/S". since dominant heartbeats 

0 can have more than one unique morphology Block 542 increments Num_Template to one to indicate one stored tem- 
plate Flow of control moves to block 545. the operation of which will be discussed later. 

Once at least one template is created block 516 is answered negatively, and block 520 slides the beat to classify 
across the first stationary template As previously discussed, the first template usually corresponds to a first heartbeat 
classification of D/S. meaning the template for dominant and supraventricular waveforms. As block 520 slides this beat 

> across the first template the minimum area difference between the beat to classify and the first template is calculated 
3lock 525 asks »f this minimum area difference is less than a threshold, if it is. the beat matches the first (D/S) template, 
and block 530 averages the new beat data into ihe template it matched. In the preferred embodiment, a weighted 
Htverage is used where the existing template is given more weight in Ihe averaging process lhan the new beat. Bk^ck 
530 also keeps track of the number of itmes a beat matched this template, as well as when a beat most recently 

' matched this template 

1 1 block 525 !S answered negatively, block 535 increments the template counter Block 538 verifies that the maximum 
number of templates to check has not been exceeded, which would indicate ail templates have been checked If block 
535 "s answered affirmatively flow of control loops back lo blocks 51 8 and 520. where the beat is slid across Ihe second 
■emplate The second template exemplarily corresponds to a classification of ventricular ectopic ("V"). Block 520 again 
computes the minimum area difference, and block 525 again asks if this minimum area difference is less than a thresh- 
ed If \{ »s. the beat matches the second (V) lemplate. and block 530 averages the new beat data into the template it 
matched if no match was found the subroutine loops through blocks 535, 536. 51 S. 520. and 525 until a match is 
Jound or uniil block 533 is answered negatively, indicating thai all existing templates have been checked for matches 
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Once a rnatch is found and block 530 averages ihe new beat data inlo the template it matched, block 537 verihes 
the classification ot ihe template As new beats are averaged into the existing templates, it is possible lor the classiti- 
caiion oi Ihe template to cnange For example a template originally classifred as "V may be reclassilied to 'D" as 
more beats are averaged into u, BIock 545 checks to see if the beat matched a D/S template. If so. an additional inquiry 
^Dcut the beat must be made before (he beat can be classified This is done in block 550 which asks if the beat was 
early II so the beat^ is classified as supraventricular ("S") in block 555, and flow of control loops back to block 501 to 
'ook to see I there are more beats to classify If not. the beat is classified as dominant ("D") in block 558 and flow of 
control looDS back to olock 501 to look to see if there are nrK)re beats to classify. If block 545 determined that the beat 
matched a template other than the D/S template block 550 classifies the beat as the classification corresponding to 
fo the template it matched For example, if the beat matched the ventricular ectopic ("V") template, the beat would be 
classified as ventricular ectopic Flow of control loops back to block 501 . as previously discussed. 

Referring again to block 536. »f block 53G is answered negatively all existing templates have been checked and 
none of tnem match Block 539 then asks ii Num jernplaie is less than Maxjernpiate - a counter inoicating ihe max- 
imum number of templates If not block 540 creates a new template for this beat. 

Block 542 classifies the new template and beat and increments Numjempiates, as drscussed previously 

If block 539 is answered negatively, the maximum number of templates has been reached Block 541 overwrites 
the template with the least recent update Lske block 542. block 543 classifies the new template, but does not increment 
Mum jemplates since the number of templates did not change Those skilled in the an will appreciate that the number 
of lemplates actually created can vary depending on the amount of noise in the environment and whether any ectopic 
^0 beats are detected, 

When block 501 determines that there are no more beats to classify, block 590 displays the classified heartbeats 
on display 66 of computing unit 60 {Fig. 2) One such exemplary display is shown in Fig. 15. The subroutine returns 
-n block 599 to block 205 of Fig. 5 

V. Beat Alignment 

Fig 9 shows a flowchart of the operation of subroutine 600 performed by alignment logic 75 of computing unit 60 
of the preferred embodiment of the invention Block 501 determines which 3 ECG waveforms are the least noisy In 
:he preferred embodiment, this is done by using the information obtained in block 301 of the Delect QflS Complexes 

JO in ECG Waveform subroutine 300. although this determination could be made directly in this subroutine by continuously 
calculating the signal to noise ratio or other indication of signal quality on the twelve ECG waveforms and ranking these 
waveforms from highest to lowest, or by using some other technique. Block 603 sets a lead counter to look at the first 
of the three leasl noisy ECG waveforms Block 605 receives ECG dala from the ECG waveform determined by the 
iead counter After the ECG data is low pass fillered in block 608. block 610 determines the dominanl template amongst 

^•5 me beats to be used to construct the representative heartbeat. In the preferred embodiment, this is done by looking 
at the number of times the templates used by the classification logic were matched by a beat, as determined in block 
530 in Fjg 3 By definition, this will be a D'S template 

Block 620 creates and stores an alignment template in storage 69 of computing unit 60 In the preferred embodi- 
ment, the alignment template is created with a oeal which matched the dominant template determined in block 61 0. 

-0 ' where a por,ion of the dominant template around the QRS is normalized . An alternate embodiment has been contem- 
plated wnere 5iep 61 0 is sKtpped. and block 620 creates an alignment template by finding the first beat classified as 
D .via beat classification logic 74 or a conveniicnai method of beat classifjcation) and using this beat as the alignment 
template 

Block 530 sets a beat counter to 1 Block 640 gets the next dominant beat for this ECG waveform and normalizes 
-5 ^ portion of this beat arounc the QRS complex. For the purposes of this application, beats classified as "O' are referred 
to herein as "dominant" -vhile beats classified -hs 'V '5^ or "Q" are referred to as "non-dominanl* In the preferred 
embodiment non -dominant beats are excludea from alignment and from the determination of a representative heart- 
beat, since {hese beats can contaminate the representative heartbeat. 

Block 650 slides this beat across the staijona/y alignment template, while computing the value of the sum of the 
absolute values of the difference, between the beat and the alignment template. This value is referred to herein as the 
area difference The position where the area difference is the minimum ( f^/linimum Area Difference) is the position where 
the beat is oesi aligned with the alignment template and this position is saved for this beat in storage 69 of computing 
unit 60 Block 6E " • :remenis the beat counter 3lcck 660 checks to see if the beat counter is greater than or equal to 
[he number of b. • to align for this ECG waveform If .-lol flow of control loops back to block 640 to gel the next 
dominant beaf If so block 670 increments the .'ead counter 

Block 675 cnecKS to see if the lead counter is qroater than or equal to 3 the number of least noisy ECG waveforms 
selected m block 601 If block 675 .s answered noqativeiy flow of conlrcl loops back to block 605 where the alignment 
process is repeated fo' the next ECG wavefcrn^ sf block 675 is answered affirmatively, block 630 corrects and stores 
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each beat lime on each ECG waveform with the median value of the atignmenl limes just determined independently 
on the three least noisy leads. This is done to minimize the effect of noise which results in beats on different ECG 
waveforms best aligning at slightly different times (i e jitter). The subroutine returns in block 699 to block 700 in Fig. 
5 Fig t6 shows a new beat being slid across a stationary template 

5 

VI. Repfeaentatlve^Beat Creation 

Fig 10 shows a flowchart of the operation of subroutine 700. performed by representative heartbeat creation logic 
77 of computing unit 60 of the preferred embodiment of the invention. Block 701 determines which beats to use in 

'0 creating the representative heartbeat. In the preferred embodiment, only "dominant" beats, aligned via tfie beat align- 
ment sleps discussed above are used. If more "dominant" beats are available than the number needed to construct 
a representative heartbeat, those with the most similar morphology are used. For example, if two or more Dominant 
iempiates were created in classification subroutine 500. only the beats that match the dominant template that contains 
the most beats wiH preferably be used.. Block 7'05 sets a counter which keeps track of which ECG waveform the 

>^ representative heartbeat is being created for to one. 

Block 710 checks to see if the lead number counter is greater than or equal to the maximum number of ECG 
waveforms If not. block 71 5 reads the aligned "dominant" beats for this ECG waveform. Bkx;k 720 sets a lime pointer 
to zero 

Block 725 gets a lime slice of data for each aligned beat at the instance of time Identified by the time pointer In 
^0 ihe preferred embodiment, this data is the magnitude of each of the aligned dominant heartbeats al this moment ol 
time. Block 730 sorts the magnitudes for this slice of time from smallest to largest. Block 735 discards a percentage 

ol Jhe smallest and largest magnitudes In a noisy environment, the beat classification subroutine may incorrectly 
classify beats as being dominant. These errois result in misciassified beats being incorrectly included in the aligned 
beats. The inmmed average done here results in a robust and accurate representative heartbeat even in the presence 
of misciassified beats and high noise present on dominant beats. In the preferred embodiment. 33% of the smallest 
magnitudes and 33°'oOf the largest magnitudes are discarded, although other'values could be used. Block 740 averages 
the remaining magnitudes for this time slice Block 750 stores the average magnitude for this time slice in a represent- 
Mive beat array in storage 69 of computing unit 60 Block 755 increments the time pointer to the next slice of lime, and 
block 750 checks to see if Ihe time pointer has reached its maximum value. If not. flow of control loops back to block 

^0 725 to determine the average magnitude of the other slices of time to complete the representative heartbeat array. Fig. 
17 shows exemplary aligned beats being time sliced using the process described above. 

When block 760 is answered affirmatively, block 765 forward and reverse filters the representative heartbeat stored 
•n the representative heartbeat array, and stores the result back into the array in block 770. An allemate embodiment 
has been contemplated where this step is skipped. Block 775 increments the ECG waveform counter and flow of 

J5 control loops back to block 71 0 to create a representative heartbeat for each of the other ECG waveforms. Once block 
710 determines thai representative heartbeats have been created and stored for each of the ECG waveforms, block 
750 displays the representative heartbeats on display 66 of computing unit 60. An exemplary display of a representative 
heartbeat is shown in Fig. 13 The subroutine retums in block 799 lo block SOO of Fig. 5. 

-0 VII. Measurements 

Fig 1 1 shows a flowchart of the operation ol subroutine BOO. performed by measurements logic 78 of computing 
unit 60 of the preferred embodiment of the invention Block 801 obtains the representative heartbeats for each of the 
ECG waveforms. In the preferred embodimeni. this is done by reading the representative heartbeat array stored m 
step 770 of Fig 10 Allemaiively. representative heartbeats created using a different method, including conventionally 
known methods for creating representative heartbeats could be used. 

Block 503 measures the earliest QRS onset and the latest QRS offset ol all of the representative heartbeats ob- 
tained in blocK 301 These values are used for many of the measurement-^^ 'hat will be made for these representative 
heartbeats, as will soon be discussed. Block 505 sets a counter that keeu . ack of the ECG waveform for which the 
representative heartbeat is being measured BlocK 5 1 0 gets the representative heartbeat for this ECG waveform. Block 
515 determines the isoelectric ievel of the representative heartbeat. - In the preferred embodiment, this is the average 
level of the 13 msec of data prior to tne earliest ORS onset. Block 820 determines the R wave amplitude for this 
representative heartbeat In the preferred embodiment ihis is the maximum positive value between the earliest QRS 
crsei and the latest QRS offset with an adjustment made to correct for elevated ST segments al the latest QRS offset 
•( :ne "T" wave is so large ihat it impacts ihe aelermiration of Ihe R wave amplitude. 

Block 525 determines the ST level in the pfeferred emoodimenl this is the average of 10 msec around the user- 
determined ST measuremeni point of the representanvo heaabeat. Block 830 determines Ihe ST slope. In the preferred 
e.'Tibodiment. ihis is determined by using a best line fit between the latest QRS offset and the ST measuremeni point 
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cf the representative heartbeat Block 535 determines Ihe ST integral. In ihe preferred embcxjiment. this is delermrr'ied 
by computing the sum of the negative area between the latest QRS offset and the ST measurement point of the rep- 
resentative heartbeat 

Block =50 updates a measurement confidence flag for each measurement taken. In the preferred embodiment, 
historical information and physiologic limits are used to set these measurement flags to either a "low" or "high" confi- 
dence These confidence flags can be displayed to a user in a variety of ways, including the term "low" or "high- 
displayed next to a irieasurement changing the color of the measurement on the display (e.g.. green means high, red 
means lowi etc A "low" confidence flag would indicate to the cardiologist or other medical professional that a meas- 
urement IS not physiologic or has changed in a non-physiologic manner and should be manually reviewed for correct- 
ness Block 555 increments the ECG waveform counter Block 560 checks to see if the ECG waveform counter exceeds 
the maximum number of ECG waveforms if not flow of control loops back to block 31 0 to repeat the process for the 
other ECG waveforms. (I so. block 360 displays the measurements on display 66 of computing unit 60 (Fig. 2). An 
exbinpiary display of these measurements, displayed alorig ^r^itli the representative heartbeat, is shown In Tig. 19. A 
cardiologist looking at the representative heartbeat and the measurements shown In Fig. 1 9 would see that there is a 
depressed ST segment, indicating that the patient undergoing a stress test has coronary artery disease. The subroutine 
returns in block 599 to block 2l0 of Fig. 5. 



Claiins 

^0 

1. A medical device ( 10) for classifying a plurality of heartbeats in an ECG waveform, said ECG waveform obtained 
from a plurality of ECG electrodes, said medical device compnsing: 

a computing unit (60) having a processor (65) memory (66) and storage (59), said computing unit further 
comprising heartiseai classification logic(74) said heartbeat classification logic further comprising: 

means for comparing a first heartbeat from said pluraiily of heartbeats with a first template stored in memory 
or storage said first template corresponding to a first heartbeat classlticalion: 

■^^ means lor determining if said first heartbeat matches said first template, responsive to said comparing means 

.f said first heartbeat matches said first template. 

means for averaging said first heartbeat into said first template. 

35 

means for classifying said first heartbeat as said first heartbeat classification ; and 
means lor stonng said first heartbeat classification for subsequent diagnosis. 
-^0 2. The meoical device of claim 1 further comprising: 

a display for displaying an indicia of said first heartbeat classification 

3. The medical device of claim i wherein said averaging means averages In a manner that is weighted in favor of 

-5 said first template 

4. The medical device of claim l said heartbeat classification logic further compnsing: 

if said first heartbeat does not match said rirst template means for comparing said first heartbeat with a second 
'^^ template stored in memory or storage corresponding lo a second heartbeat classification: 

means for determining if said lirsi heartbeat matches said second template, responsive to said comparing 
means, and 

;t said first heartbeat matches said second template means for classifying said first heartbeat as said second 
heartbeat classification 

5. The medical device of claim 4 said heartbeat clnssification logic further compnsing. 
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means for averaging said first heartbeat into sard second template. 

6. The medical device of claim 5. wherein said averaging means averages in a manner that is weighted in favor of 

Scii6 second template 

7. The meaicai device of claim 4 said heanbeat classification logic further comprising: 

if said first heartbeat does not match said second template, means for creating a third template for said first 

heartbeat and 

means for classifying said third template as a third heartbeat classification. 

8. The medical device of claim 1 wherein said first heartbeat matches said first classification and said first classifi- 
cation IS a dominant/supraventricular ectopic classificaiiorv said medical device further comprising: 

menns for checking (o see if said first heanbeat was early: 

it said checking means determined that said first heartbeat was early, means for classifying said first heartbeat 
as supraventricular ectopic, and 

if said checking means determined lhat said first heartbeat was not early, means for classifying said first heart- 
beat as dominant 

9. The medical device of claim 4 said heartbeat classification logic further comprising: 

if said first heartbeat does not match said first template or said second template, means for determining whether 
said first template or said second template was last updated the longest time ago: 

means for creating a third template for said first heartbeat: 

means for overwriting said secono template with said third template, responsivd to said determining means 
determining that said second template was updated the longest time ago: and 

means tor classifying said third template as a third heartbeat classification. 

10. The medical device of claim 2. wherein said first heartbeat matches said first classification and an indicia of a 
ventricular ectopic classification is displayed on said display. 

11. The medical device of claim 2. wheiein said first heartbeat matches said first classification and an indicia of a 
questionable classification is displayed on said display 
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